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Key Objectives 


To study the 29-30 May 2012 deep convective storm observed 
during the Deep Convective Clouds and Chemistry (DC3) 
experiment over Oklahoma, including its: 

- Convective transport of trace gases 

- Associated lightning occurrence and nitrogen oxide (NO x ) production 

Simulate the observed storm using WRF-Chem 

Compare the physical features of the simulated storm against 
aircraft and ground-based observations 

Add flash rate parameterization schemes (FRPSs) to the model and 
identify the best match to observations 

Determine NO production scenario for 1C and CG flashes following 
lightning-generated NO x (LNO x ) scheme 


Background 


Storm system developed ~21Z May 29 
along KS/OK border and continued until 04Z 
May 30 

Aircraft sampled storm and its environment 
from 20Z May 29 to 01Z May 30 

— DC-8 focused on storm inflow and outflow 

— GV and Falcon concentrated on outflow 



N EXRAD Composite Reflectivity 2240Z on 29 May 


Ground-based instrumentation included: 

- Dual-Doppler radar (NEXRAD level II 
regional) 

- National Lightning Detection Network 
(NLDN) cloud-to-ground flash data 

- Oklahoma Lightning Mapping Array (LMA) 
flash initiation density data 



Blue circles: LMA stations 

Green outline: Extent of 3-D lightning mapping capability 
Gray outline: Extent of 2-D lightning detection 
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WRF-Chem Model V3.5 


• Nested domains: 15-km and 3-km 

• Initialized with DART and GFS for 
boundary conditions 

• Used coarsely prescribed IC:CG 
ratios (Boccippio et a!., 2001) 


LOCATION OF WRF PROCESSORS 
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4- Grid Cell Center LONGITUDE (°W) □ Processor 


Type of Scheme 

Selection for Simulation 

Microphysics 

Morrison 

Planetary boundary layer 

Mellor-Yamada-Janjic (MYJ) 

Radiation 

Rapid radiative transfer model for GCMs (RRTMG) 

Flash rate 

Maximum vertical velocity (W max ) 

Lightning-generated NO x (LNO x ) 

DeCaria et al. (2000, 2005) 








Flash Rate Parameterization Schemes 


• Based on simulated thunderstorm's physical features 

• Six types have previously been used in cloud- 
resolving models: 


Type of FRPS 

Equation (flashes min' 1 ) 

References 

Maximum vertical velocity 

5.0 X 10-6 x Wmax 4.5 

Price & Rind, 1992 

Cloud top height 

3.44 x 10 5 x H 4 9 

Price & Rind, 1992 

Updraft volume 

6.75 xl0- n xw 5 - 13.9 

Deierling & Petersen, 2008 

Ice water path 

33.33 x|WP- 0.17 

Petersen et al., 2005 

Ice mass flux product 

9.0 x 10 15 x (f p x f np ) + 13.4 

Deierling, 2006; 
Deierling et al., 2008 

Precipitation ice mass 

3.4 x 10 8 x p m - 18.1 

Deierling et al., 2008 






LNO x Parameterization Scheme 
(DeCaria et al., 2005) 

Gaussian vertical distribution 
of 1C (bimodal) and CG (single 
mode) NO production based 
on typical lightning flash 
channel distributions 

Lightning channels set to 
maximize at -15°C (CG and 1C) 
and -45 °C (1C) 

500 moles NO per 1C and CG 
flash (Ott et o\., 2010) 

Horizontal placement of NO 
based on reflectivity > 20 dBZ 
in each grid cell 




Methodology 


• Compared flash rate trends over the 
observed (OK LMA and NLDN) and 
model-simulated storm's lifetime 

• Created moving spatial masks at 
10-min intervals for comparison of 
observed and model-simulated storms 

• Used offline calculations, with 
adjustment factors, to analyze the 
six FRPS trends 

• Calculated NLDN total flashes given 
NLDN CG flashes and mean IC:CG ratio 
for the storm region (3.9 ± 0.49), which 
is based on Boccippio et al. ( 2001 ) 



GOES-13 1-km Visible at 0008Z (NCAR/EOL) 
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Comparison of Storm Features from 

Initial Simulation 


Model-Simulated Composite Reflectivity 2300 Z on 30 May 



Black rectangle represents the spatial mask surrounding the 
model-simulated cell of interest at 2300Z May 30 


NEXRAD Composite Reflectivity 0010Z on 30 May 



D 13 X 4* CO 73 


Reflectivity (dBZ) 

Red rectangle represents the spatial mask surrounding the cell of 
interest at 0010Z May 30 


Model-simulated storm: 

• Began ~1-1.5 hour before observed storm 

• Exceeds area of observed storm by roughly a factor of two 
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Compared instantaneous 
flash rates from WRF at 
10-min intervals with 
corresponding 1-min 
periods from the observed 
NLDN flash rates 

Model-simulated flash 
rate trends are adjusted 
90 minutes later to 
coincide with observations 


Adj. instant 10-min WRF output vs. instant 10-min Adj. NLDN Obs. 
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Observations 

WRF - precipitation ice mass 
WRF - updraft volume 
WRF - maximum vertical velocity 
WRF - ice mass flux product 
WRF - ice water path 
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Flash Rate 

Parameterization Scheme 

Total Flashes 
Prior to Scaling 

Scaling Factor 

Maximum vertical velocity 

3,951 

1.1310 

Cloud top height 

708 

6.3138 

Updraft volume 

21,118 

0.2116 

Ice water path 

4,452 

1.0035 

Ice mass flux product 

36,745,336 

0.0001 

Precipitation ice mass 

164,749 

0.0271 


Flux product, precipitation ice mass, and 
updraft volume trends are similar to the 
increasing trend of observations 

Timing of W max and ice water path peaks is 
similar to observations (140, 200, & 310 min) 

Magnitude of observed primary peak greater 
than those in FRPSs 

W max and ice water path schemes need the 
least adjustment to match observed total 
flashes at each 10-min interval (4,468 
flashes) 






LMA Flash initiation density 

DDHHMM = 300120 ODHHMM = 300130 
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Count per grid pixel per 10 min 


10 min total 
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W max FRPS overestimat es the total flashes of both the NLDN (~2.5) and LMA (~3.5) 
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* Plots courtesy of M. Belo 


Observations taken in cloud-free air to the 
south of the storm system 

To define the real background air at 9-12 
km, the 10 th percentile is used to remove 
any influence from old convective outflow 


Aircraft 

Time of UT Background 
Sampling (Z) 

DC-8 

20:40-21:10 (black) 

GV 

22:15-22:30 (red); 22:58- 
23:05 (yellow) 


Obs. 

6-9 km (DC-8) 

9-12 km (DC-8, GV) 

Species 

Median (ppb) 

10 th percentile (ppb) 

NO 

0.025 

0.029 

no 2 

0.034 

0.033 

N0 X 

0.059 

0.062 
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Model-simulated vertical cross-section 
taken in cloud-free air to the south of the 
storm system 


Aircraft 

Time of UT Background 
Sampling (Z) 

Model 
Time (Z) 

DC-8 

20:40-21:10 

20:00 

GV 

22:15-22:30; 

22:58-23:05 

21:00; 

22:00 


Altitude 

NO 

(PPb) 

no 2 

(PPb) 

NO x 

(PPb) 

6-9 km 

Obs. 

0.025 

0.034 

0.059 

(Median) 

Model 

0.020 

0.017 

0.01-0.05 

9-12 km 

Obs. 

0.029 

0.033 

0.062 

(10 th percentile) 

Model 

0.041 

0.023 

0.05-0.5 


Model-simulated NO x at 9-12 km is given 
as a range. Although the observed NO x is 
found at the lower end of this range, it 


should be kept in mind that the observed 
value represents the 10 th percentile. 
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Initial LNO x Analysis from Aircraft 


NEXRAD Composite Reflectivity 2230Z on 29 May 



Cloud-free 
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Radar observations overlaid with 10-min intervals 
of aircraft NO x measurements 

NO x measurements increase from sampling in 
cloud-free upper tropospheric air (0.02-0.1 ppbv; 
upper left) to making transects through anvil 
outflow (peaks to ~2 ppbv; bottom left) 

Model-simulated total LNO x a magnitude lower in 
anvil outflow than observations (bottom right) 
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2300Z on 29 May 
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New Model Simulation 


• Changes include: 

— One domain (1 km) vs. two-way nested domain (15/3 km) 

- NAM 18Z vs. DART/GFS 12Z initialization 
— CLM vs. NOAH land surface scheme 

• Model-simulated storm now begins about an hour after 
observed storm 

• Area of observed and model-simulated storm is roughly 
similar 

• Comparison of offline WRF flash rate calculations against 
observations indicates an overestimation of a factor of ~10 
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Conclusions 

• Based on offline calculations, W max FRPS was selected for use in model: 

- Needs little adjustment to match the observed total flashes 

- Coincides with several of the observed flash rate peaks 

• Scale up model-simulated flash rates in offline calculations and scale down online: 

- May partly be due to how offline calculations are computed 

• Model overestimate of observed flashes may be due to: 

- Area of model-simulated storm ~2x larger than observed 

- Observed storm passes over northern edge of LMA 

• Initial look at NO x chemistry in UT air undisturbed by storm: 

- At 6-9 km, NO values similar between aircraft and model, and model-simulated N02 
underestimates observations by ~0.02 ppbv 

- At 9-12 km, the 10 th percentile NO x values are similar (~0.06 ppbv) between the aircraft and 
model 

• New simulation at finer resolution produces model-simulated storm of roughly 
same size as observed 

• WRF-Chem model estimates of P(NO x ) in works 

— Model-simulated LNO x underestimates aircraft observations in anvil outflow by roughly a 
magnitude 

— LNO x calculations are being evaluated in the model and in offline calculations for verification 
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Future Work 


• Perform a trace gas simulation 
and analysis of NO x , CO, and 
0 3 using WRF-Chem 

• Compare model-simulated 
LNO x against aircraft measured 
NO x 



• Determine NO production 
scenario per 1C and CG flash 
that best matches aircraft 
observed NO x mixing ratios 



Investigate 0 3 changes 
downwind of flight 
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Photo by C. Cantrell 




NEXRAD Composite Reflectivity 0000Z on 30 May 
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Model-simulated total LNO x 
about a magnitude lower in 
anvil outflow than aircraft 
observations 


Model-Simulated LNOx Composite Reflectivity 2300Z on 29 May 


CG LNO x (pptv) at 7 km 
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1C LNO x (pptv) at 11 km 


Total LNO x (pptv) at 11 km 
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NEXRAD Composite Reflectivity 0000Z on 30 May 






Model-simulated total LNO x 
about a magnitude lower in 
anvil outflow than aircraft 
observations 


Model-Simulated LNOx Composite Reflectivity 2300Z on 29 May 


CG LNOx (pptv) at 7 km 
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1C LNOx (pptv) at 11 km 


Total LNOx (pptv) at 11 km 
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